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Abstract 


Porous wall wind tunnels have been used for several decades and have proven 
effective in reducing wall interference effects in both low speed and transonic testing. 
They allow for testing through Mach 1, reduce blockage effects and reduce shock wave 
reflections in the test section. Their usefulness in developing computational fluid 
dynamics (CFD) codes has been limited, however, by the difficulties associated with 
modelling the effect of a porous wall in CFD codes. Previous approaches to modelling 
porous wall effects have depended either upon a simplified linear boundary condition, 
which has proven inadequate, or upon detailed measurements of the normal velocity near 
the wall, which require extensive wind tunnel time. 

The current work was initiated in an effort to find a simple, accurate method of 
modelling a porous wall boundary condition in CFD codes. The development of such a 
method would allow data from porous wall wind tunnels to be used more readily in 
validating CFD codes. This would be beneficial when transonic validations are desired, 
or when large models are used to achieve high Reynolds numbers in testing. 

A computational and experimental study was undertaken to investigate a new 
method of modelling solid and porous wall boundary conditions in CFD codes. The 
method utilized experimental measurements at the walls to develop a flow field solution 
based on the method of singularities. This flow field solution was then imposed as a 
pressure boundary condition in a CFD simulation of the internal flow field. The 
effectiveness of this method in describing the effect of porosity changes on the wall was 
investigated. Also, the effectiveness of this method when only sparse experimental 
measurements were available has been investigated. The current work demonstrated this 
approach for low speed flows and compared the results with experimental data obtained 
from a heavily instrumented variable porosity test section. 
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The approach developed was simple, computationally inexpensive, and did not 
require extensive or intrusive measurements of the boundary conditions during the wind 
tunnel test. It may be applied to both solid and porous wall wind tunnel tesjs. 
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Chapter 1 


A 


Introduction 


1.1 Background 

Ventilated wall wind tunnels have been in use for several decades and have been 
useful in reducing wall interference effects at subsonic and transonic speeds and 
allowing for testing through Mach 1. A series of improvements have been made to the 
earliest ventilated wind tunnels leading to modem porous wall test sections. The 
usefulness of porous wall wind tunnels for computational fluid dynamics (CFD) 
validation efforts has been limited, however, by difficulties associated with modelling 
the porous wall boundary condition. 

It has long been recognized that the corrections to wind tunnel data for open and 
closed test sections were of opposite signs [1]. Furthermore, in transonic testing, shock 
waves which impinge on a solid boundary are reflected as shock waves, whereas shocks 
which impinge on a free air boundary in an open jet test section are reflected as 
expansion waves. Theodorsen suggested that a wind tunnel might be constructed which 
would reduce wind tunnel wall interference by using a partially open boundary condition 
in which one, two or three walls would be removed from a solid wall test section [2]. 
Wright and Ward tested one of the first successful ventilated wall wind tunnels which 
used several stream wise slots in the tunnel walls [3]. They found that blockage 
interference was reduced with this wind tunnel. In addition, they found that ventilated 
walls alleviated choking problems at transonic speeds and permitted testing through 
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Mach 1. These two effects have led to ventilated wall test sections being used for both 
low speed and transonic work [4]. The work by Wright and Ward led to other ventilated 
wind tunnels, also using streamwise slots to ventilate the test section walls. It was soon 
realized, however, that streamwise slots allowed for the reflection of shock and 
expansion waves from the walls. These reflected waves could impinge on the model and 
result in data which was uncorrectable [5]. 

Porous wall tunnels alleviated this problem by significantly reducing the shock 
reflections from the walls [6]. The porous wall was constructed with a pattern of small, 
discrete holes in the wall. The differential resistance wall was a further refinement to the 
porous wall concept, using holes with axes inclined to the normal. The differential 
resistance wall was found to provide cancellation of both shock and expansion waves 
and allowed for inflow and outflow resistance of the wall to be balanced [1]. 

1.2 Motivation for Modelling the Porous Wall Boundary Condition 

As high speed computers have increased the flexibility and capability of 
computational modelling of flow fields, a new emphasis has been placed on obtaining 
wind tunnel data which may be used to calibrate and validate computational fluid 
dynamics (CFD) codes [7]. Increasingly, there has been a trend toward modelling entire 
wind tunnel flow fields, including support struts and wind tunnel walls [8, 9]. This has 
led to an increased use of solid wall wind tunnels in transonic testing. The simplicity of 
modelling a solid wall boundary condition has made their use attractive despite the 
disadvantages of substantial wall interference [10]. 

Porous wall wind tunnels have long been used for both low speed and transonic 
wind tunnel testing because of their desirable effect of reducing wind tunnel wall 
interference. With mounting evidence of the importance of Reynolds number on flap 
gap sizing and maximum lift predictions [11, 12], the need for performing high 
Reynolds number tests to validate CFD codes which will be used in the design process is 
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of growing importance. While a solid wall wind tunnel may be used for such tests, wall 
interference can become a limiting factor on the size of models. 

In general, for test sections of identical dimensions, a porous wall wind tunnel 
will allow testing of larger models at transonic speeds without shock reflections from the 
tunnel walls impinging on the model. In addition, porous walls may eliminate shock 
boundary layer interactions on the walls and reduce such interactions on the model by 
eliminating reflected shocks from the walls. Thus, the use of porous wall data in CFD 
validation may reduce the grid refinement required near the wind tunnel walls while 
providing a more realistic assessment of the applicability of a turbulence model to free 
air calculations. Ultimately, the goal of many CFD validation efforts is to develop CFD 
codes which are capable of predicting free air performance. These codes may be useful 
for design purposes even if they are not sophisticated enough to resolve multiple shock 
boundary layer interactions which may be associated with shock reflections from solid 
wall wind tunnels. 

Some attempts have been made to model discrete slots in CFD codes either by 
modelling the slots in an approximate manner or by solving coupled equations to 
describe the boundary condition [13, 14]. However, porous wall wind tunnels have been 
shown to be superior to slotted walls at cancelling shock wave reflections [6]. Porous 
walls pose a problem in CFD validations, however, since the modelling of the wall 
boundary is more difficult for a porous wall wind tunnel. 

The intricacies of the porous wall make it prohibitive to depict the actual wall 
geometry in a viscous CFD grid. Since the porous wall wind tunnel may have several 
thousand holes on the walls, modelling the individual holes and the viscous effects 
associated with each hole in a CFD grid is not possible given the current limitations on 
computer speed and memory. Thus, the effect of the porous walls must be dealt with 
either by correcting the test data to free air conditions, or by modelling the porous wall 
by appropriate means in the CFD code. 
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1.3 Previous Approaches to Modelling Porous Walls 

Numerous approaches for correcting porous wall data to free air conditions have 
been proposed. These methods use a variety of approaches, based on model pressure 
and force measurements [15], wall boundary pressure or velocity measurements, or 
pressure rail measurements [4, 16]. These methods generally produce a global 
correction to the velocity and angle of attack, based on classical reflection techniques 
[1]. The test data is then taken to be equivalent to data from a test in free air at the 
corrected angle of attack and velocity. Additional corrections to drag and moment 
coefficients and Mach number are sometimes included [1]. This results in a useful 
comparison of bulk flow measurements, such as lift and drag coefficients. However, 
these methods are of limited use in performing CFD validations and calibrations because 
they have the effect of altering the entire flow field. When performing CFD validations, 
it is desirable to compare as much of the flow field as possible [7]. 

Other approaches have sought to develop a boundary condition which may be 
used at the porous wall boundary in the CFD code. These methods have depended on 
either a simple, universal boundary condition, [1,6] or on detailed measurements of the 
boundary conditions during the wind tunnel test which are then imposed as boundary 
conditions in the CFD code [17-21]. 

A universal boundary condition has proven difficult to determine experimentally, 
and to implement computationally. The porosity parameter contained in such a 
boundary condition can only be determined by extensive dedicated testing. In addition, 
this parameter may be different for otherwise identical top and bottom walls, and it may 
depend on Mach number, stagnation pressure, model size and orientation [22]. 
Additionally, such a boundary condition may be destabilizing in CFD codes and can 
actually prevent convergence [23, 24]. 

Most current efforts in ventilated wall interference research have been directed 
toward making detailed measurements of velocity or pressure in the wind tunnel and 
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using these to develop a boundary condition in CFD codes. The approach used by King 
and Johnson in modelling a slotted wall boundary condition was to use a rail suspended 
midway between the model and the tunnel wall to make pressure measurements at every 
point in the flow corresponding to a boundary point in the CFD code [17]. Later, to 
allow for additional grid refinement, a spline fit was made to this data as a means of 
interpolating pressure measurements [18]. This approach required extensive 
measurements of the pressure on planes away from the wall and has not been extended 
to three dimensions. 

Jacocks used detailed measurements of static pressure, normal velocity adjacent 
to the wall and boundary layer thickness at the wall to develop an equivalent inviscid 
normal velocity profile based on boundary layer integration along the wall [19]. Crites 
and Rueger extended this method for limited three-dimensional problems by using 
extensive pressure measurements on the walls and extensive calibrations of the tunnel 
walls to estimate the normal velocity through the porous wall and, again, develop an 
equivalent inviscid normal velocity profile. This profile was then used as a boundary 
condition in numerical solvers for the purpose of developing incremental corrections to 
wind tunnel data [20, 21]. These approaches required detailed measurements of the 
wall boundary conditions and large amounts of dedicated test time for calibration at each 
plenum pressure, Mach number and Reynolds number, making them costly in practice. 
In most practical cases, the extensive measurements and calibrations required at each 
test condition have limited these techniques to two-dimensional flows. 

1.4 Current Approach 

The approach presented here represents an attempt to describe the effects of a 
porous wall in a CFD code based on sparse measurements of the flow field in the wind 
tunnel. The usefulness of this theory has been demonstrated by carrying out wind tunnel 
experiments in a heavily instrumented, low speed (M s .07) research wind tunnel and 
performing CFD simulations of the experiments. This test section utilized variable 
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porosity walls and allowed for extensive pressure and velocity measurements to be made 
near the boundaries of the test section. 

The theoretical approach to describing the porous wall boundary utilized the 
classical porous wall boundary condition as a means of interpolating and extrapolating 
boundary conditions measured during wind tunnel tests. A solution for the flow near the 
walls, based on the method of singularities, was developed using pressure measurements 
made on the walls during wind tunnel tests. The singularity expressions used in this 
method satisfied the classical porous wall boundary condition. The method of 
singularities solution was then used to specify the pressure at CFD boundary point 
locations. 

The current work extended the classical method of singularities by including 
higher order singularities. The current work also differs from classical approaches in 
that it attempted to model the boundary condition in the CFD code, rather than applying 
corrections to the data based on measurements of the boundary condition. Furthermore, 
unlike some limited previous attempts at modelling the boundary condition, the current 
work utilized experimental pressure measurements on the walls and used classical 
porous wall theory to extend the measurements. This approach did not require extensive 
calibrations of the wind tunnel walls, and used static pressure measurements on the 
walls, which may be obtained rapidly during wind tunnel tests. 

In two dimensions, this approach was developed using analytic expressions for 
the singularity strengths. In three dimensions, a numerical solution was required to 
determine the pressure profiles induced by singularities in the presence of the porous 
wall boundary condition. Once obtained, however, these numerical solutions may be 
used for multiple tests. 

The current work consisted primarily of two-dimensional, low speed studies. 
Experimental work was conducted in the Stanford Low Speed Wind Tunnel and 
computational studies were performed using an incompressible flow solver. Although 
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the Mach number range of the research tunnel used for this work was very low, some 
consideration was given to the implementation of this method for higher speed flows. 
Limited consideration was also given to the extension of this method to three- 
dimensional flows. 
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Chapter 2 


Porous Wall Theory 


2.1 Introduction 

One of the major advantages of porous wall wind tunnels is their ability to reduce 
or eliminate shock wave reflections in transonic testing. In addition, they have the effect 
of reducing wall interference for low speed flows, making them useful for cases where 
substantial wall interference is expected in low speed testing [25]. Classical approaches 
to correcting data from porous wall wind tunnels are slowly being replaced by efforts 
aimed at modelling the porous wall boundary condition in CFD codes. 

2.2 Development of Classical Porous Wall Boundary Condition in Two Dimensions 

It has widely been assumed that the flow through the porous walls is basically 
viscous in nature. This has led to a simple theory, based on an analogy to pipe flow, for 
describing the porous wall boundary condition [1,6]. This boundary condition assumes 
that the airfoil perturbations are sufficiently small near the walls so that linearization 
applies up to high subsonic Mach numbers near the walls. An implicit restriction of this 
theory is that the flow at the walls must be subsonic. In practice, however, flow near the 
walls has been found to be subcritical at freestream Mach numbers up to .9, lift 
coefficients as high as 1.5 and a height to chord ratio of 4 [26]. 

Classical porous wall theory [1] typically uses a point source, point vortex, and 
point source and vortex doublets to describe the free-air flow past a lifting airfoil, as 
given in Equation 2.1. 
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Note that Equations 2.4 and 2.5 may be interpreted as first derivatives with respect to x 
of Equations 2.2 and 2.3, respectively. 

The potential <j>p is a solution of the linearized potential equation: 
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An interference potential, which accounts for the interference caused by the 

presence of the walls must be solved for. Again, in the region between the walls, 
^ < y < (jty must also satisfy the linearized potential equation: 
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The free air potential and interference potential together form the potential 
function <|>: 
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which is required to satisfy the following boundary conditions: 


P 


u 


+ $1 - o 

dx dy 


at y 


h 

2 


Pl 


d§ 

dx 


3<t> 

dy 


= 0 at y = — 


( 2 . 10 ) 

* 
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The basis for these boundary conditions is easily seen. The x-axis is aligned with the 

freestream flow in the wind tunnel so that = v, and^ = u. According to linearized 

ay ox 

small perturbation theory, u should be proportional to the difference between free stream 
and local pressure. If the plenum is maintained at freestream static pressure, this will 
also be the pressure difference across the wall. Based on an analogy to pipe flow, the 
velocity normal to the wall, v, should be proportional to the pressure difference existing 
across the wall. The constant of proportionality is the porosity parameter, Py for the 
upper wall and Pl for the lower wall. Different porosity parameters may be used for the 
upper and lower walls since the flow through the two walls is different in character. In 
the case of the upper wall, there is a suction over most of the wall resulting in flow into 
the test section. This flow comes from the plenum where the flow is essentially at rest. 
For the lower wall, there is an increased pressure, resulting in flow out of the test 
section. This flow has considerable streamwise momentum. The result is that flow 
through the lower wall often separates within the holes in the wall [1, 20]. Thus, 
identical upper and lower walls may offer differing resistance to flow through the walls. 
In part, this may be compensated for by inclining the holes in the wall at an angle with 
respect to the normal direction. The different porosity parameters specified on each wall 
offer another means of accounting for this difference. 

Analytic expressions for the velocities induced by simple singularities in the 
presence of the boundary conditions of Equations 2.10 and 2.1 1 have been developed [1, 
27] and are summarized here. 
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The wind tunnel coordinates are nondimensionalized as. 


X =7T- Y = r 

|5h h 

(2.12) 

and the complex position function Z is defined as. 


Z = X + i Y 

(2.13) 

and the complex velocity notation is defined as. 
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W = Ww + W F = H - 

(2.16) 


With the form of the free-air potential due to a singularity known, it remains to 
determine the form of the complex velocity function W^r. 

The complex velocity funtion may be found by the method of images. The 
approach of Reference 1, which is followed here, determines the interference potential 
arising from a singularity placed between two infinite boundaries described by 
Equations 2.10 and 2.11. As this approach may be found in detail in the literature [1, 
27], only the result will be shown here. 

Considering the potential flow due to singularities in the presence of an infinite 
porous wall boundary condition as described by Equations 2.10 and 2.11, it may be 
shown that the complex velocity functions for a source, W a , and a vortex, Wy, are: 


W c = | [(B + B) + X(Pu)X(Pl) ] 

(2.17) 

W y = iy(B-E) 

(2.18) 
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where: 


*iu±ik (z -z 0 ) 

_ 1 e L 
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0< Hf + *L <1 (2.19)' 


E = - — 
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0 < iU±iL < x (2.20) 


X(Pu)X(P L ) = 


1 P„*P t =0 
0 Otherwise 


( 2 . 21 ) 
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( 2 . 22 ) 


Velocity Functions for source and vortex doublets in the presence of infinite porous 
walls may be obtained by differentiation of the above formulas: 


W„ = — (B + E) 

P P h dZ 


W (B = ^ 4 (B-E) 


P h dZ 


(2.23) 

(2.24) 


Higher derivatives of the source and vortex may also be taken if additional singularities 
are desired. In this work, up to the fourth derivatives of the source and vortex were 
retained in the singularity expansions. 

It is noteworthy that the velocities induced by a source in this formulation are 
discontinuous at P=0 (see Equations 2.17 and 2.21). Classical wind tunnel interference 
corrections model the wall interference using infinite strip theory. This approach 
develops corrections based on linear theory solutions for simple singularities in the 
presence of two infinite walls. When these walls are solid, the mass introduced by the 
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source must exit either upstream (x = or downstream (x = +°°.) Generally, a bias is 
introduced into the velocity functions in order to force the velocity disturbance to go to 
zero at x = -co (see Equation 2.21), but the influence of the source far downstream 
remains non-zero. However, for a porous wall, no matter how small the porosity, the 
mass contribution from the source may leak out through the walls. Thus, for an infinite 
porous wall, the velocity contribution from a source will go to zero at both x = and 
x = +oo. This discontinuity is well known [28]. The use of infinite strip theory 
introduces finite test section errors [29], but it is still common in practice as an 
engineering approximation. The effects of the discontinuity at P=0 caused by the 
infinite strip theory will be considered further in Chapter 5. The velocity profile on the 
wind tunnel inflow plane and the presence of a breather section at the outflow plane also 
introduce finite test section effects. In this work, singularity solutions were developed 
using infinite strip theory. Finite test section effects associated with the inflow and 
outflow plane were minimized by specifying boundary conditions based on experimental 
measurements on these planes. A velocity profile, based on five-hole probe 
measurements was specified on the inflow plane and a constant pressure was specified at 
the breather location. 

2.3 Methods of Describing the Effect of a Porous Wall in CFD Codes 

Classical porous wall theory has been widely used to apply corrections to the 
measured wind tunnel data based on the interference potential <t>w The use °f sucil 
corrections is undesirable for CFD validations, however, because they alter the entire 
flow field in order to impose a correction which is valid at only one point. Recent 
efforts in CFD validations have attempted to model wall boundary conditions directly in 
the computations. For a solid wall wind tunnel test, this may be accomplished by 
specifying either a slip (inviscid) boundary condition at the edge of the boundary layer, 
or a no-slip (viscous) boundary condition at the location of the walls. For a porous wall. 
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however, the complexity of the wall geometry does not allow for the wall to be modelled 
exactly given the current limitations on computer memory and speed. 

Some attempts have been made to use the boundary conditions of Equations 2.10 
and 2. 1 1 directly in CFD calculations. This approach is questionable for three reasons. 
First, there is not a good theoretical approach for determining the value of the porosity 
parameter contained in these equations and so the porosity parameter must be 
determined from experimental data [1]. Second, there is experimental evidence that the 
linear relationship between perturbation velocities implied by Equations 2.10 and 2.11 
do not apply over the entire wall [19]. Third, and most important, the application of 
Equations 2.10 and 2.11 as boundary conditions in CFD codes is destabilizing and can 
actually prevent convergence [23, 24]. Still, classical porous wall theory does provide 
some insight into the nature of the flow near a porous wall. 

Some attempts have also been made to model a porous wall boundary condition 
by specifying a normal velocity profile on the walls. The normal velocity profile 
developed in these approaches is an equivalent inviscid normal velocity profile based on 
an integral boundary layer analysis. This analysis requires a knowledge of the normal 
velocity at the wall surface. These approaches have been successful in describing the 
effects of a porous wall, but their application has been limited by the extreme 
requirements for wind tunnel calibrations and measurements during tests. These 
approaches require either detailed normal velocity measurements near the wall during a 
test or detailed pressure measurements on the wall coupled with extensive calibrations of 
the wind tunnel walls in order to determine the relationship between pressure differences 
across the wall and velocities induced through the wall [19, 20, 21]. The calibrations 
must be performed at every Mach number and Reynolds number which will be used 
during testing and they require an adjustable pressure plenum. Either implementation of 
this approach requires extensive wind tunnel tests for the purpose of determining 
boundary conditions. Furthermore, the integral methods used depend upon very detailed 
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measurements during the test, limiting the facilities in which they may be used. While 
some attempts have been made to apply these procedures to simple three-dimensional 
cases, their usefulness is limited by test facility requirements. , 

2.4 Theoretical Approach of the Present Work 

The approach used for modelling the porous wall boundary condition in this 
research was to make use of classical porous wall theory as a means of interpolating and 
extrapolating limited experimental measurements of the pressure on the porous walls. 

In this work, the classical porous wall boundary condition of Equations 2.10 and 
2.1 1 was used to develop pressure profiles on the walls using the method of singularities. 
The singularity solutions were developed as discussed in Section 2.2. Singularity 
strengths were determined based on least squares matching of velocity profiles inferred 
from pressure measurements made on the walls. The resulting method of singularities 
solution was then used to specify a pressure boundary condition in the CFD code. Thus, 
flow field corrections were eliminated in preference to a means of specifying 
experimental boundary conditions in the CFD simulation. Furthermore, by using 
singularity solutions based on the classical porous wall boundary condition, 
experimental data could be interpolated and extrapolated over the length of the walls. 

Classical two-dimensional porous wall theory typically uses four singularities (a 
source, a vortex, a source doublet and a vortex doublet) to develop corrections to the 
flow field. The current approach extended classical theory by using as many as ten 
singularities. These singularities included a source, a vortex, and first, second, third and 
fourth derivatives of the source and vortex in the streamwise direction. Furthermore, in 
the current apporach, no corrections were applied to the wind tunnel data. Instead, the 
pressure profile existing on the walls during the experiments was modelled in the CFD 
code. In addition, classical porous wall boundary conditions were used only to develop 
the singularity solutions, and were not applied directly in the CFD code. In this work. 
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the porosity parameter was found by minimizing the root-mean-square errors in 
matching pressure coefficient measurements on the walls. 

2.5 Extension to Three Dimensions 

Extension of this approach to three dimensions requires a panel method, or other 
numerical solver, to determine three-dimensional singularity solutions in the presence of 
a porous wall boundary condition. Again, as in the two-dimensional case, simple 
singularities and their derivatives may be used. In this work, a line source of span 2yQ, 


<t>S=-^-ln 


P(y-yo)+\ 

f x 2 +p 2 (y-y 0 ) 2 + (pz) 2 

P(y+yo) + -\| 

' x 2 + p 2 (y + yo ) 2 + (P z ) 2 


and horseshoe vortex of span 2yQ, 


(2.25) 


<t>v = 


4k 


A tan 




4 - A tan 


;(y+yo) 


z(x 2 + P 2 (y + yo ) 2 + P 2 z 2 f 2 


- Atan 


x(y-yo) 


z(x 2 +p 2 (y-y 0 ) 2 +p 2 z 2 j 


yM 


(2.26) 


along with the first, second and third derivatives with respect to x of the line source and 
horseshoe vortex, the first, second and third derivatives with respect to yo of the line 
source and horseshoe vortex, and a first derivative with respect to z of the horseshoe 
vortex were used to develop the method of singularities solution. A uniform freestream 
term was also included in each case as a means of biasing the source term. In Equations 
2.25 and 2.26 above, x is the Cartesian coordinate in the streamwise direction, y is the 
Cartesian coordinate in the spanwise direction, and z is the Cartesian coordinate in the 
vertical direction. The singularities are located at x = 0, z = 0, and centered about y = 0. 
They have a span of 2yo- 
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The panel method used constant strength doublet panels to model the wind tunnel 
walls [30]. The walls were modelled over 24 chord lengths. Geometric stretching of the 
panels in the streamwise direction was used to allow for better resolution of the solution 
near the location of the model. A total of 384 panels were used to model each of the 
four walls. Inflow and outflow planes were not panelled. This approach, which models 
only a finite portion of the walls rather than treating them as infinite, results in 
continuous behaviour at P = 0. Figure 1 shows the arrangement of panels on the walls. 
In Figure 1, as in Equations 2.25 and 2.26, the singularities are located at x = 0, z = 0 
and centered at y = 0. The span of the singularities was equal to the wing span. A 
simple shift in the coordinate system may be used to place the singularities at any 
streamwise location desired. For the calculation of singularity strengths based on least 
squares matches to experimental data, the singularities were shifted to the quarter-chord 
location on the wing. 

For three-dimensional testing, the boundary conditions analogous to Equations 
2.10 and 2.11 are: 


3<f> 3<(> 

dx dz 

= 0 

at 

h 

z = — 
2 

(2.27) 

3<|> 3<j> 

dx dz 

= 0 

at 

h 

z = 

2 

(2.28) 

9<j> + 3<t> 
dx dy 

= 0 

at 

y=W 

(2.29) 


In the case of full-span testing, an additional boundary condition of the form of Equation 
2.29 may be imposed on the remaining wall. In the case of half-span testing, as in this 
work, the remaining side wall is a plane of symmetry. The boundary conditions of 
Equations 2.27-2.29 were enforced at control points located at the center of each panel. 
Bicubic spline fits were used to interpolate solutions between control points on the 
panels. 
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Chapter 3 


Two-Dimensional Governing Equations and Navier 
Stokes Solver 


3.1 Introduction 

For numerical simulation of the wind tunnel tests performed in the Stanford Low 
Speed Wind Tunnel, an incompressible Navier-Stokes solver (INS2D) was used [31]. 
Although methods of describing porous walls in CFD codes are of primary interest for 
transonic codes, an incompressible code was used in this case because of the very low 
Mach number of the tests conducted in this research tunnel. The method described in 
this and the previous chapter may also be used with transonic flow solvers [32]. 

3.2 Governing Equations 

The equations presented here have been nondimensionalized as follows: 

~ _ u i ~ _ x i z _ t 11 ref ~ _ P ~ Pref 

u i — x i ~ L ~ P _ 2 

u ref x ref x ref P u ref 

ty — — — ^ — v = — - — = Re _1 (3.1) 

P u ref x ref u ref 

where the ~ superscript denotes nondimensional quantities. Hereafter, the ~ superscript 
will be dropped for convenience. 

In the above equations, ui = u or v, the Cartesian velocity components for i = 1, 
2, respectively and xi = x or y, the Cartesian spatial coordinates for i = 1, 2, 
respectively. The pressure, density and kinematic viscosity are denoted by p, p, and v. 
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respectively. The subscript re f denotes reference quantities; uref is the measured 
freestream velocity upstream of the model, x re f is taken as the airfoil chord. 

The Navier-Stokes equations for incompressible, two-dimensional, constant 
density flow may then be written in conservative form as follows: 


3u dv 
dx dy 


= 0 


9u 

di 






where: 


(3.2) 
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(3.5) 



(3.6) 


For turbulent flows, these equations represent Reynolds averaged quantities. The 
Boussinesq approximation is used for the Reynolds stress. The viscous stress tensor 
may then be written as: 


' c ij=( v+v t) 


du* duj 
3xj 3xj 


(3.7) 


where v t is the turbulent eddy viscosity. 

As the thrust of this work is the development of a methodology for describing 
porous wall boundary conditions in CFD codes, the details of the CFD code and 
turbulence modelling will not be discussed here. A brief description of the CFD code is 
included in the following sections. Additional details about the CFD code may be seen 
in References 31 and 33. The results shown from the INS2D code used a Baldwin-Barth 
turbulence model. The Baldwin-Barth turbulence model is in wide use in a variety of 
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CFD codes currently. Additional details about this turbulence model may be seen in 
Reference 34. 

3.3 Numerical Solver 

The INS2D code uses a primitive variable formulation of the incompressible 
Navier-Stokes equations. In the primitive variable formulation, the equations are written 

with pressure and velocity as the dependent variables. 

The equations are transformed into generalized coordinates using the 

transformations: 

% = $(x, y, t) 

T|=r|(x,y,t) ( 3 - 8 ) 

The system of equations then becomes: 



(3.9) 



where J is the Jacobian of the transformation and: 


(3.10) 
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The contravariant velocity components, U and V are given by: 

U = £ x u + £ y v V = q x u + T|yV (3.13) 

The metrics of the transformation are represented with subscript notation as, for 


example: 



The INS2D code uses an artificial compressibility approach for solving the 
incompressible Navier- Stokes Equations [31, 35]. In this approach, a time derivative of 
pressure is added to the continuity equation (Equation 3.9) and an artificial 
compressibility factor, P, is introduced: 

^ = - p(V'0) (3-14) 

dr 

This results in a hyperbolic system of equations. This system may then be marched to a 
steady-state solution in pseudo-time. The right hand side of Equation 3. 10 and the 
divergence of the velocity field approach zero as the solution approaches steady-state. 

This hyperbolic system of equations also has artificial pressure waves which are 
finite in speed. Therefore, compressible flow algorithms may be used for solving the 
system of equations. The INS2D code uses third order accurate, locally upwind 
differencing of the convective fluxes and second order accurate central differencing of 
the viscous terms. The upwind differencing is biased by the signs of the eigenvalues of 
the local flux Jacobian, following the method of Roe [36]. The set of numerical 
equations is solved using a Gauss-Seidel type line relaxation scheme. 

All computations were performed on the CRAY Y-MP computer at NASA Ames 
Research Center. Solutions typically converged in approximately 1500 iterations, 
requiring approximately 30 minutes of CPU time on the CRAY Y-MP computer. 
Convergence histories for cases using slip wall boundary conditions were nearly 
identical to those for cases using a specified pressure boundary condition. 

Additional details concerning this algorithm and its numerical implementation 
may be found in Reference 31. 

3.4 Boundary Conditions 

The method of artificial compressibility used in the INS2D code allows for 
boundary conditions to be specified using the method of characteristics, similar to 
compressible flow solvers. In the case of the artificial compressibility scheme, however, 
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the method of characteristics is merely a mathematical construct which allows for 
waves to exit at the boundaries of the computational domain without reflection. 

There are three characteristic lines associated with three eigenvalues of the 
system of equations: u, u+c, u-c, where c is the artificial speed of sound in this method. 
At an inflow boundary, there is one characteristic line carrying information from the 
interior of the solution to the boundary. At an outflow boundary there are two 
characteristic lines carrying information from the interior of the solution to the 
boundary. 

The approach used in this work was to treat both the outflow plane and top and 
bottom walls as outflow boundaries. Although not strictly correct for a wall boundary 
which allows both inflow and outflow, this approach is consistent with that used by 
other researchers [17-21]. The problems associated with making velocity measurements 
near a wind tunnel wall usually result in only one flow variable being specified as a wall 
boundary condition in wind tunnel simulations. Additionally, an abrupt change from 
inflow to outflow boundary conditions on the wall usually results in discontinuities in 
the solution. 

The outflow boundary condition, as implemented in the INS2D code, can most 
easily be understood by consideration of a limiting case in which the method is applied 
to a uniform, orthogonal grid. In this case, for an outflow boundary, the method reduces 
to: 
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where the n-subscript denotes differentiation in the direction normal to the boundary. 
For other cases, the metrics of the transformation and local rotations also affect the 

equations. 

With p specified at the boundaries, and for cases where the artificial 
compressibility parameter P » u, v, as was the case in this work, this method results in 
an extrapolation of interior velocities to the boundary. It does, however, still allow for 
small variations in pressure and velocity normal to the surface to avoid wave reflections. 

It therefore allows for faster convergence than a simple extrapolation scheme. 

3.5 Two-Dimensional Grid Generation 

Grids were developed which conformed to boundaries on which experimental 
data was available. The inflow plane of the grid corresponded to the inflow plane on 
which five-hole probe data was taken. Top and bottom boundaries of the grid 
corresponded to the nominal (0-degree inclination) top and bottom wall locations. The 
outflow plane corresponded to the breather section location. These grids were used for 
all slip wall computations and for all cases in which pressure was specified as a 
boundary condition. A special grid was developed for the viscous solid wall case. For 
that case, the top and bottom boundaries of the grid corresponded to the inclined wall 
locations. 

All grids used with the INS2D code were developed using GRIDGEN software 
[37]. Surface definitions were based on spline fits to approximately 2000 measurements 
of the model geometry. The very high resolution of measured model geometry allowed 
for excellent geometry definition in the grids as well as resolution of model defects 
associated with the attachment of the model flap and spoiler. Orthogonality was 
maintained at the airfoil surface by using exponential blending of grid lines near the 
surface. The Thompson Middlecoff elliptical solver was used within the GRIDGEN 
software to develop the grids. This solver maintains boundary distributions of grid 
points in the interior of the grid and allows for easy refinement of the grid in boundary 
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layers over the airfoil. Spacing was controlled so that approximately 15 to 20 grid points 
were located within the boundary layer on the airfoil. All of the grids developed were 
C-grids. Normal grid dimensions were 250 x 70, although a test case was run using a 
350 x 105 grid to examine the effect of grid refinement. An example of these grids may 
be seen in Figure 2. 



Figure 2. Grid used for two-dimensional CFD simulations for a = 5 degrees. 


25 


Chapter 4 


6 

Experimental Apparatus 


4.1 Wind Tunnel Test Section 

A specially constructed test section was built for use in the Stanford Low Speed 
Wind Tunnel. This test section was heavily instrumented with static pressure tappings 
on the walls and with five hole probe access on the inflow and outflow planes and on 
planes near the top and bottom walls. The test section was designed to allow for two- 
dimensional testing and three-dimensional half-span testing. For two-dimensional 
testing the top and bottom walls were porous and both side walls were solid. For three- 
dimensional half span testing, the model was mounted to a solid wall and the remaining 
three walls were porous. 

The test section was .457 m x .457 m and the porous portion of the test section 
was .76 m long. The porous walls were constructed using a layered design to allow the 
open area ratio to be varied. The layer of the wall nearest the flow was a slotted layer 
with .00254 m streamwise slots located at .0254 m intervals. The middle layer of the 
wall consisted of inserts which had .0229 m by .00254 m rectangular holes located at 
.0254 m intervals. These holes were inclined to the normal at 15 degrees, giving the 
porous wall a differential resistance to inflow and outflow. The outer layer of the wall 
was also slotted similar to the inner layer and could be slid laterally to vary the open area 
ratio of the wall. Open area ratios of the wails could be varied from a solid wall 
condition to an open area ratio of 9-percent. An atmospheric plenum was used for all 
porous wall tests. An exploded view of this wall design may be seen in Figure 3. 
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An atmospheric breather section was also used in all tests and was located at the 
end of the test section. The nominal frees tre am velocity for all tests was 24. m/sec and 
the nominal Mach number for all tests was M = .07. The nominal Reynolds number 
based on chord was Re = 3.2 x 10^ for the airfoil tests and 2.4 x 10^ for the three- 
dimensional wing tests. The test section top and bottom walls were adjustable and were 
inclined normal to the freestream to compensate for boundary displacement layer growth 
through the test section. Displacement layer thickness was measured by a boundary 
layer rake located near the end of the porous portion of the test section. Displacement 
layer thicknesses were measured on both the top and bottom walls with the model in 
place. The rake was removed during normal tunnel operation. A schematic of the wind 
tunnel may be seen in Figure 4, and a photograph of the test section may be seen in 
Figure 5. 

4.2 Wind Tunnel Models 

The two-dimensional model used in these tests was a Boeing Advanced 
Transport Airfoil. A cross section of this airfoil may be seen in Figure 6. The model 
had a flap and a spoiler which were in the retracted position for these tests. The airfoil 
chord was .203 m and the maximum thickness ratio was 11.3 percent Measurements of 
the model were made with a Leitz Precision Measuring Machine with an accuracy of 
2.54 x 10~6 m, or approximately .00125 percent of the airfoil chord. Approximately 
2000 measurements were made to define the airfoil cross section. Grit was applied to 
the airfoil between 2.5 and 5.0 percent of the chord [38, 39]. Nominal grit diameter was 
.03175 cm. Oil flow studies were undertaken to ensure the two-dimensionality of the 
flow field and to ensure that the grit applied to the model successfully triggered 
transition on the airfoil. 

The three-dimensional model used in these tests was an unswept half-span wing 
with an aspect ratio of 1.5, a span of .229 m and a chord of .152 m. The wing had a 
constant NACA 0018 airfoil cross section. The wing tip was a half body of revolution 
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with diameter equal to the local thickness of the wing. Grit was applied to the wing 
between 2.5 and 5.0 percent of the chord. Nominal grit diameter was .03175 cm. Oil 
flow studies were undertaken to ensure that the grit successfully triggered transition on 

a 

the wing. 

4.3 Test Section Instrumentation 

The test section was instrumented with static pressure tappings on the porous 
walls and five-hole probe access near the top and bottom walls and on the inflow and 
outflow planes. The tunnel was instrumented with 78 static pressure tappings on the 
centerlines of the upper and lower walls and a total of 482 static pressure tappings on 
the upper and lower walls. An additional 121 pressure tappings were located on the 
porous side wall. Pressure tappings were also located on both models. The five-hole 
probe was used to determine flow inclinations and velocity perturbations on the inflow 
and outflow planes as well as planes near the top and bottom walls. 

Static pressure data was acquired by three Scanivalves which were calibrated on 
each run. Together, these Scanivalves were capable of making 576 pressure 
mesurements on each run. Nominal freestream velocity was determined by static 
pressure measurements taken in the inlet contraction bell. These measurements were 
sufficiently far upstream so that interference from the model was negligible. The 
velocity measurements were calibrated in empty tunnel runs. 

4.4 Wind Tunnel Flow Quality 

Empty tunnel flow quality surveys were conducted and may be seen for solid 
wall and porous wall cases in Figures 7 and 8, respectively. A velocity variation of 
approximately 4 percent may be seen in these surveys at the inflow plane location. This 
effect appears to be local to the inflow plane of the test section, and is much less 
pronounced near the model location. Instrumentation access holes restricted the range of 
positions at which measurements could be made. In these cases, as in the data which is 
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shown later, the nominal free stream velocity is based on static pressure measurements 
made in the contraction section of the wind tunnel. 

4.5 Error Analysis 

All pressure measurements were obtained with instrumentation calibrated against 
a secondary standard traceable to the National Bureau of Standards. Pressure 
measurements on the walls and airfoil were made with three Scanivalves, each having 
four 48-port barrels. Two ports on each barrel were used for a reference and calibration 
pressure. This allowed the pressure transducer in each barrel to be calibrated on each 
run. A settling time of approximately 1 second was allowed after each advancement of 
the Scanivalve barrel. A low pass (10 Hz) filter was used to filter data from the pressure 
transducers. Pressure readings were based on an average of 200 measurements taken at 
1 kHz. Standard deviations were computed for these measurements and were typically 
0.01-percent or less of the mean data values. A calibration pressure of 0.100 psi (.689 x 
103 N/m 2 ) was used for wall pressure measurements and a calibration pressure of 0.200 
psi (1.378 x 10 3 N/m 2 ) was used for airfoil pressure measurements. Uncertainties in the 
calibration pressure measurements were ± 0.6 x 10*3 psi (4.13 N/m 2 ), based on 
manufacturer supplied data. Additionally, discretization by the data acquisition board 
resulted in a minimum resolution error of .1 x 10*3 psi (.69 N/m 2 ) for the airfoil 
measurements and .05 x 10*3 p S i (34 N/m 2 ) for the wall measurements. Non-linearity 
of the pressure transducers was estimated at less than . 2 -percent based on a comparison 
of known pressure differences and transducer outputs. These errors combined to give an 
overall error of approximately .9-percent in Cp measurements on the airfoil and 1.2- 
percent in Cp measurements on the walls. 

Five-hole probe pressure data was obtained by dedicated 1 psi Druck pressure 
transducers. These transducers were calibrated on each run using the secondary standard 
discussed above. Calculation of flow angularity without direct measurements of velocity 
perturbations led to results of unknown accuracy. The relationship between pressure 
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coefficient measurements with a five-hole probe and flow angles is complex. Adding to 
uncertainties in pressure mesaurements were uncertainties regarding probe alignment. A 
qualitative assessment of the probe data could be made by comparing data taken from 
identical runs, or by observing the variations in repeated probe measurements made at 
the same point in the flow field. This qualitative assessment suggested that typical probe 
error was on the order of 0.5 degree. 
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Figure 3. Exploded view of variable porosity wall. 
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Figure 4. Schematic diagram of Stanford Low Speed Wind Tunnel. 



Figure 5. Variable porosity test section. Airfoil model is shown installed. 
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Figure 6. Boeing Advanced Transport airfoil cross section. 



Figure 7. Empty tunnel, solid wall flow quality surveys at test section entrance and 
model location. 
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Figure 8. Empty tunnel, porous wall flow quality surveys at test section entrance and 
model location. 
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Chapter 5 


A 


Results 

5.1 Two-Dimensional Case 

Both computational and experimental studies were undertaken to evaluate the 
method of singularities for application to two-dimensional flows. Preliminary 
computational studies used a CFD solution to determine the number of measurements 
and the number of singularities required to develop an adequate description of the flow 
near the walls using the method of singularities. The experimental data was used to 
determine the suitability of this method to actual flow fields. Effects considered were: 
1) the number of singularities and flow field measurements required to describe the flow 
near the walls, 2) the ability of this method to capture the effect of changes in the open 
area ratio of the walls, 3) the effects of uncertainties in the porosity parameter specified 
in the method of singularities, 4) the effects of singularity placement, 5) comparisons of 
alternate methods of describing a solid wall boundary, 6) the effects of grid refinement 
in the CFD code, 7) the effects of using sparse data to develop a wall pressure profile, 
and 8) flow inclinations near the boundary, as computed by the CFD code and measured 
by the five-hole probe. 

5.1.1 Preliminary Computational Studies 

As a means of evaluating the method of singularities, a CFD simulation using 
slip wall boundary conditions was developed. This simulation was then treated as a 
psuedo-experiment and data from the wall boundaries was used to develop a method of 
singularities solution. The use of a computational result instead of experimental data for 
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this purpose allowed for a detailed comparison of pressure profiles on the walls and a 
comparison of the effects of using sparse data to develop the method of singularities 
solution. Figure 9 shows the pressure profiles on the walls as predicted by the CFD , 
simulation using slip wall boundary conditions for an angle of attack of 5 degrees. 

Figure 10 shows the root-mean-square errors resulting when various numbers of 
singularities and data points were used to develop the potential flow solutions. The data 
points were distributed as uniformly as grid spacing would allow over the top and 
bottom walls. Up to four derivatives of a source and vortex were retained in the method 
of singularities solutions. The root-mean-square errors shown in Figure 10 were based 
on a point-by-point comparison of the pressure predicted by the CFD simulation and the 
pressure predicted by the method of singularities at each grid point on the upper and 
lower walls. The root-mean-square errors have been normalized by the peak pressure 
coefficient on the walls. As may be seen in Figure 10, a relatively small number of 
singularities was adequate to describe the far field flow. In addition, a small number of 
measurement locations was sufficient to develop a potential flow solution. Additional 
data from the walls allowed for a better description of the pressure profile in the method 
of singularities, but beyond 25 data points the improvement was very small. 

The use of additional singularities resulted in the need for additional data points 
in order to produce acceptable fits to the velocity profiles. It was found that the number 
of data points required to produce a good fit to the velocity profiles was approximately 
twice the number of singularities retained in the potential flow solution. 

While these trends were observed for a variety of airfoil simulations [32, 40, 41], 
caution should be exercised when applying this method to different geometries. Some 
computational experimentation, of the type just described, may be required in order to 
develop confidence in this method for different geometries such as multi-element 
airfoils, blunt bodies, or configurations with very low height-to-chord ratios. 
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5.1.2 CFD Simulations of Experiments 

The next step in evaluating this method was to use it to perform simulations of a 
series of wind tunnel tests. The tests were conducted in the variable porosity test section 
described in Chapter 4. The pressure profiles measured on the walls during the 5 degree 
angle of attack tests may be seen in Figure 11a. This figure shows data from both the 
solid wall tests and the porous wall tests with a 9-percent open area ratio. Also shown in 
Figure 11a are the potential flow solutions based on this data. These potential flow 
solutions were developed using the 77 measurements on the wind tunnel walls and a 
total of 10 singularities (a source, a vortex, and the first four derivatives in the 
streamwise direction of a source and vortex.) Figures 12a-17a show similar results for 
angles of attack of 6, 7, 8, 9, 10, and 11 degrees. As can be seen in these figures, the 
potential flow solutions showed excellent agreement with the experimental data. 

The potential flow solutions shown in Figures 11 a- 17a were imposed as 
boundary conditions in the Navier-Stokes solver and the resulting pressure profiles 
obtained on the airfoil are shown in Figures llb-17b. Also shown in Figures llb-17b 
are the airfoil pressure profiles obtained in the wind tunnel tests. Experimental and 
computational results compared very well. The largest discrepancies were 
approximately 6-percent of the peak values. These discrepancies occurred at the 
location of the grit transition strip on the model. At locations away from the grit strip, 
however, agreement was much better with typical discrepancies of approximately 1- 
percent of peak values. A comparison of airfoil pressure profiles from the solid wall and 
porous wall tests showed that this method simulated the effect of a porous wall very 
well. As the porosity of the wall was varied, the trend in pressure profiles on the airfoil 
was duplicated in the CFD simulations by imposing the potential flow solution as a 
pressure boundary condition in the CFD solver. 

Also of interest in Figures llb-17b are the slight deviations seen in the pressure 
profiles on the airfoil between x/c = 0.75 and x/c = 0.80. These deviations may be seen 
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on both the upper and lower surfaces and are present in both the experimental and CFD 
data. They are the result of small ridges on the airfoil model between the main airfoil 
and the retracted flap and spoiler. 

The root-mean-square errors in matching experimental pressure coefficients on 
the walls using the method of singularities are summarized in Table 1. Also shown in 
Table 1 are the root-mean-square errors between measured pressures on the airfoil and 
pressures predicted by the CFD code with the method of singularities boundary 
condition imposed. In both cases, these root-mean-square errors are based on a point- 
by-point comparison at all locations where experimental data was available. The root- 
mean-square errors shown for wall and airfoil cases have been normalized by peak 
pressure coeffiecients measured on the wall and airfoil, respectively, in each test. 
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Table 1. Root-mean-square errors in matching pressure coefficients on the walls using 
the method of singularities, and in matching experimental pressure coefficients on the 
airfoil with the INS2D CFD code. 

5.1.3 Effect of Changes in Porosity Parameter 

In the results shown above, the porosity parameter used in the method of 

singularities was determined by minimizing the root-mean-square errors in matching the 

wall pressure profile. As a means of evaluating the sensitivity of this method to 

uncertainties in the porosity parameter specified in the method of singularities, a variety 

of porosity parameters were used to develop potential flow solutions in the method of 
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singularities. Figure 18 shows root-mean-square errors in matching pressure coefficients 
on the walls as a function of the porosity parameter specified in the method of 
singularities for the 9-percent porous wall, 8-degree angle of attack case. The root- 
mean-square errors shown in this figure have been normalized by the peak pressure 
coefficient measured on the walls. These results are typical of the results at other angles 
of attack. It is noticeable that the curves shown are discontinuous at P = 0. The reasons 
for this discontinuity have already been discussed in Chapter 2. In all the cases 
examined, the optimum porosity parameter was found to be P = 0. For non-zero 
porosities, the optimum porosity depends on the number of singularities retained and the 
optimum tends to be very shallow. This again points out the danger of using Equations 
2.10 and 2.1 1 directly as a boundary condition in CFD codes. 

Figures 19 and 20 show a comparison of results obtained by specifying a 
porosity parameter of 0. and 1. in the method of singularities solution. These results are 
for the porous wall tests at an angle of attack of 8 degrees, and using 10 singularities in 
the method of singularities. Figure 19 shows the resulting pressure coefficient profiles 
on the walls, and Figure 20 shows the resulting pressure coefficient profiles predicted by 
the CFD code on the airfoil. As can be seen in Figure 19, the finite test section effects 
associated with the change in porosity are most noticeable near the inflow and outflow 
planes of the test section. The trend in pressure profiles was captured with either 
specified porosity. As can be seen in Figure 20, the discrepancies in the pressure 
profiles of Figure 19 have little effect on the airfoil pressure profiles. Overall, this 
method is relatively insensitive to uncertainties in the porosity parameter. By using a 
least-squares approach to determine singularity strengths, this method allows singularity 
strengths to vary as needed in order to provide an optimum fit to the experimental data 
for any specified value of porosity. This effect may be seen in Figure 21, which shows 
the singularity strengths in the method of singularities as a function of the porosity 
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parameter specified on the walls. For clarity, only 6 singularity strengths are shown in 
this figure. 

5.1.4 Effect of Changes in Singularity Placement 

In the results discussed above, all singularities were placed at the quarter-chord, 
tunnel centerline location (x/c = 0.25, y = 0.0). Figure 22 shows the effect of varying the 
location of the singularities. Singularity locations were varied from x/c = 0.0 to x/c = 
1.0. The resulting root-mean-square errors in matching the pressure profile on the walls 
for the solid wall, 8-degree angle of attack case are shown. Again, the root-mean-square 
errors have been normalized by the peak pressure coefficient on the wall. As can be 
seen, this method is relatively insensitive to the position of the singularities, provided 
that the singularities are placed ahead of the mid-chord point An increase in the root- 
mean-square errors occurs as the singularities are placed further downstream, although 
even these errors are relatively small. As noted earlier, some caution should be 
exercised when applying these results. When using this method with geometries 
significantly different from those used here, such as multi-element models or extremely 
low height-to-chord ratios, it may be necessary to distribute singularities in order to 
adequately model the pressure profile on the walls. 

5.1.5 Comparison of Alternate Boundary Conditions for Solid Wall Simulations 

As a means of comparing this approach to other methods of specifying boundary 
conditions in CFD codes, CFD simulations were performed of the solid wall 8-degree 
angle of attack test using 3 alternate boundary conditions. These alternate boundary 
conditions were: 1) A slip-wall boundary condition (v = 0) applied at the nominal wall 
location, 2) A no-slip boundary condition (u = v = 0) applied at the physical wall 
location, and, 3) A free-air boundary condition (u = u©o, v = Voo, p = poo) applied at a 
distant boundary with corrections applied to the airfoil angle of attack and the freestream 
velocity. 
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The grid used for the slip wall simulation was the same grid used to specify the 
pressure boundary condition. The outer boundary of this grid corresponded to the 
nominal wall location. As discussed earlier, the walls were inclined to adjust for the 
growth of the displacement layer through the test section. 

Figure 23 shows the grid used for the no-slip wall simulation. This grid models 
the wall inclinations through the test section and uses additional refinement near the 
walls to resolve the wall boundary layer. 

Figure 24 shows the grid used in the free-air simulations. For the free-air 
boundary condition, classical incidence and velocity corrections were applied to the 
wind tunnel data [42]. The CFD simulation was performed with the airfoil at the 
corrected angle of attack and the velocity corrections were applied in non- 
dimensionalizing the data. The outer boundary was located 30 chord lengths from the 
airfoil. 

Figure 25 shows the pressure profiles on the airfoil resulting from these 
simulations, as well as the simulations using the pressure boundary condition from the 
method of singularities. Experimental data is also shown for reference. The specified 
pressure boundary condition, based on the method of singularities, produced the closest 
match to the experimental data. The slip wall and no-slip wall showed similar results 
which were also in close agreement with the experimental data. The free-air simulation 
showed significant discrepancies (approximately 18-percent of peak values) with 
experimental data at the pressure peak, although overall agreement was much closer. 
Overall agreement between these methods was close, although differences may clearly 
be seen in the area near the pressure peak on the airfoil. The overall agreement of these 
approaches demonstrated the consistency of this method with other means of specifying 
boundary conditions for solid wall tests. This agreement gave added confidence in this 
method when it was used for simulating porous wall boundary conditions. 
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5.1.6 Grid Refinement Effects 

The effects of grid refinement were also considered for the solid wall data from 
the 8-degree angle of attack case. The grids used in the above results had already been 
refined to the point where additional grid refinement appeared to have little effect on the 
solution. For comparison purposes, however, an additional result was computed on a 
further refined grid using approximately twice the number of grid points as were used in 
the standard grids typical of Figure 2. Figure 26 shows a 350 x 105 grid which was 
used to simulate the 8-degree angle of attack, solid wall case. Figure 27 shows a 
comparison of the results obtained for this case on the 350 x 105 grid and the 250 x 70 
grid which was used in the other simulations. As can be seen in this figure, differences 
between these results are minor. Some additional refinement was evidenced near the 
flap and spoiler ridges on the aft portion of the airfoil, but overall results indicated that 
the 250 x 70 grid was sufficiently refined so that remaining errors associated with grid 
refinement were minor. 

5.1.7 Effects of Using Sparse Data 

The effect of using sparse data to develop potential flow solutions was also 
investigated. Figure 28a shows 8 data points selected from the complete set of wind 
tunnel data for the solid wall test at 5 degrees angle of attack. A potential flow solution 
based on these 8 data points is also shown. In addition, the complete set of wind tunnel 
data and the potential flow solution based on the complete set of data is also shown for 
reference. For the sparse data (8 data point) solutions, only first derivatives of the 
source and vortex were retained in the method of singularities, resulting in a total of 4 
singularities being used to develop the potential flow solution. Figures 29a-34a show 
similar results for solid wall tests at 6-11 degrees angle of attack and Figures 35a-41a 
show similar results using data from the 9-percent open area ratio porous wall tests. The 
potential flow solutions based on 8 data points show close agreement with the potential 
flow solutions based on 77 data points. Root-mean-square errors between experimental 
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measurements and method of singularities solutions based on 8 data points are shown in 
Table 2. These root-mean-square errors are based on a point-by-point comparison of the 
entire set of 77 experimental measurements with the method of singularities solutipn 
based on 8 measurements. The root-mean-square errors are normalized by the peak 
pressure measured on the walls for each case. As expected, these errors were higher 
than the cases in which 77 data points were used to develop the solution (see Table 1), 
however the trend of the pressure profiles is still captured using very sparse data. This 
agreement indicates that this method was robust enough to allow for a good description 
of the boundary conditions on the walls even when only very sparse data was available 
from the wind tunnel walls. 
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Table 2. Root-mean-square errors in matching pressure coefficients on the walls using 
the method of singularities based on 8 experimental measurements, and in matching 
experimental pressure coefficients on the airfoil with the INS2D CFD code. 

The pressure profiles shown in Figures 28a-41a were imposed as boundary 
conditions in the CFD code and the resulting pressure profiles on the airfoil are shown in 
Figures 28b-41b. For reference, the airfoil pressure profiles which resulted by imposing 
the boundary conditions of Figures lla-17a in the CFD code are also shown in Figures 
28b-41b. This allows a comparison of the effect on the airfoil of using sparse data to 
develop the wall boundary conditions. As can be seen in these figures, the method of 
singularities provided a means of describing boundary conditions which was sufficiently 
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robust so that little effect was seen in the airfoil pressure coefficient profiles from 
significantly reducing the number of measurements used to develop the boundary 
conditions. Root-mean-square errors between experimental data on the airfoil and the 
CFD results using boundary conditions based on 8 data points are also summarized in 
Table 2. Again, these errors have been normalized by the peak pressure measured on the 
airfoil in each case. These errors are comparable to those for the cases in which 77 data 
points were used to develop the boundary condition (see Table 1), and indicate that this 
method may produce close agreement between computed and experimental 
measurements on the airfoil even when only very sparse data is available to develop 
boundary conditions. 

5.1.8 Flow Inclinations Near the Boundaries 

Five-hole probe measurements were made on a plane near the boundaries, as 
described in Chapter 4. These measurments were used to develop a qualitative 
assessment of the flow inclinations near the walls. 

Figures 42 and 43 show the normal velocities measured near the walls during 
porous wall tests of the airfoil at angles of attack of 5 and 8 degrees, respectively. In 
these figures, the normal velocity is plotted against the streamwise location. Velocity 
profiles are shown for both the top and bottom walls, and all velocity measurements are 
normalized by freestream velocity. Also shown in these figures are the computed 
normal velocities from the CFD code. These normal velocities are taken from a plane in 
the CFD result corresponding to the five-hole probe measurement plane. These CFD 
results were obtained by specifying the pressure boundary conditions of Figures 1 1 and 
14 in the CFD code. The difficulties of making experimental flow inclination 
measurements with the five-hole probe may be seen in Figure 43 which shows 
discrepencies between data obtained from two consecutive runs. These uncertainties 
make any quantitative comparison of computed and measured results questionable. A 
qualitative comparison of computed and measured results shows that the CFD code with 
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a pressure boundary condition is capturing the trend of the normal velocities near the 
walls over most of the boundary. 

5.2 Three-Dimensional Case 

Limited consideration was also given to three-dimensional applications of this 
method. No CFD simulations were performed for the three-dimensional data. However, 
the method of singularities was used to develop potential flow solutions based on the 
data measured on the walls. The effects of using different numbers and arrangements of 
measurements were considered, as well as the effects of using different numbers of 
singularities. 

5.2.1 Method of Singularities Comparison to Experimental Data 

Figures 44a, 44b, 45a, 45b, and 46 show the pressure profiles on the upper, lower 
and side walls from the three-dimensional solid wall tests of the half-span wing at an 
angle of attack of 20 degrees. Also shown in these figures are the pressure profiles 
developed by the method of singularities using 342 data points and 16 singularities to 

match the data. Figures 47-49 show similar results for the 9-percent open area ratio 
porous wall tests. 

As can be seen in these figures, the agreement of the method of singularities fit to 
the experimental data is better in two dimensions than in three dimensions. This is in 
part due to the finite test section effects near the inflow plane and outflow plane. The 
abrupt change in boundary conditions associated with the breather section and inlet of 
the wind tunnel test section caused flow anomalies which the method of singularities 
could not adequately model. Only data from the central region of the test section is 
shown in Figures 44-49. The root-mean-square errors based on a comparison of 
experimental pressure profiles and those predicted from the method of singularities are 
5. 1 -percent of the peak pressure coefficient for the solid wall case and 7.6-percent of the 
peak pressure coefficient for the porous wall case. Over the central region of the test 
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section, the method of singularities captured the general trend in pressure profiles on all 
three walls. 

5.2.2 Effects of Sparse Data and Reduced Numbers of Singularities 

The effects of using reduced numbers of singularities and sparse data in the three 
dimensional method of singularities were also investigated. Figures 50a-c show the 
position of the 342 pressure tappings on the walls which were used to develop the 
method of singularities solutions shown in Figures 44-49. Figures 51a-c show 223 
measurements arranged along the centerlines of each wall and three vertical or spanwise 
lines on each wall. These data points were also used to develop a method of 
singularities solution. The resulting solution was then evaluated based on a point-by- 
point comparison to the complete set of 342 experimental measurements. Figures 52a-c 
shows 136 measurements arranged along the centerlines of each wall and a single 
vertical or spanwise line on each wall. Figures 53a-c shows a total of 38 sparse data 
locations taken from the data in Figure 51a-c. The measurements indicated in Figures 
51, 52 and 53 were also used to develop method of singularities solutions which were 
evaluated based on point-by-point comparisons to the complete set of 342 
measurements. 

Figures 54 and 55 show the root-mean-square errors resulting from a point-by- 
point comparison of the entire data set and the method of singularities solutions based on 
the data at the positions shown in Figures 50-53. Results are shown in these two figures 
for the solid and porous wall, respectively. Figures 54 and 55 also show the effects of 
using various numbers of singularities to develop these solutions. The results shown are 
based on 3, 7, 11, and 16 singularities being retained in the solution. Table 3 shows the 
order in which terms were included in the singularity solutions. All root-mean-square 
errors shown in Figures 54 and 55 have been normalized by the peak pressure 
coefficients measured on the walls. 
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Number of 
Singularities 

Singularities Used 

3 Singularity Case 

Uniform freestream. Line Source, Horseshoe Vortex 

7 Singularity Case 

Uniform freestream. Line Source, Horseshoe Vortex 
First Derivatives of Source and Vortex in x-direction 
First Derivatives of Source and Vortex with respect to vn 

1 1 Singularity Case 

Uniform freestream. Line Source, Horseshoe Vortex 
First Derivatives of Source and Vortex in x-direction 
First Derivatives of Source and Vortex with respect to yo 
Second Derivatives of Source and Vortex in x-diiection 
Second Derivatives of Source and Vortex with respect to vn 

16 Singularity Case 

Uniform freestream. Line Source, Horseshoe Vortex 
First Derivatives of Source and Vortex in x-direction 
First Derivatives of Source and Vortex with respect to yo 
Second Derivatives of Source and Vortex in x-direction 
Second Derivatives of Source and Vortex with respect to yo 
Third Derivatives of Source and Vortex in x-diiection 
Third Derivatives of Source and Vortex with respect to yo 
First Derivative of Vortex in z-direction 


Table 3. Singularities included in three dimensional method of singularities solutions. 


The span of the line source, horseshoe vortex and their derivatives was equal to 
the span of the wing. Although not actually a singularity, the freestream term was 
retained in all cases to allow for the upstream contribution of the source term to be 
cancelled. As can be seen in Figures 54 and 55, the effect of adding additional 
measurements was to produce a better fit to the overall pressure profile on the walls. 
Additionally, as seen in the two-dimensional results, the addition of extra singularities in 
the potential flow solution required additional data points in order to provide an 
improved fit to the data. Overall lower root-mean-square errors were obtained for the 
cases of sparse data by using only 7 singularities instead of 11 or 16 singularities. As 
more data was used, however, the root-mean-square errors were reduced by using 
additional singularities to describe the flow field. These results suggest that 
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approximately 10 data points should be used for each singularity retained in the method 
of singularities. As with the two dimensional results, however, caution should be 
exercised when applying these trends to different geometries. Some experimentation 
with numbers of singularities and data points may be required in order to develop 
confidence in this method for geometries and test conditions which are significantly 
different from those in this research. 

Examples of the wall pressure profiles obtained by using reduced data to develop 
the method of singularities solution may be seen in Figures 56-61. Figures 56-58 show 
the pressure profiles on the top, bottom and side walls from the solid wall 20 degree 
angle of attack test and from the method of singularities using 7 singularities and 38 
measurements at the positions shown in Figure 53. Figures 59-61 show similar results 
for the 9-percent open area ratio tests. As can be seen in these figures, even when only 
sparse data was available, the method of singularities captured the general trend in the 
pressure profiles over the central region of the test section. 

5.2.3 Effects of Uncertainties in Porosity Parameter 

The effect of the porosity parameter specified in this method was insignificant. 
As discussed earlier, the three-dimensional formulation of the method of singularities is 
continuous in its behavior at P = 0. Figure 62 shows the effect on root mean square 
errors of varying the porosity parameter from 0.0 to 1.0. These root-mean-square errors 
were based on a point by point comparison of the experimental measurements from the 
20 degree angle of attack porous wall tests and the method of singularities solutions 
developed using porosities of 0.0, 0.5 and 1.0 in the panelling solutions. The complete 
set of 342 experimental measurements were used for these comparisons, and 16 
singularities (see Table 3) were used to develop the method of singularities solutions. 
All root mean square errors were normalized by peak pressure coefficients measured 
during the test. Virtually no change was observed in the root-mean-square errors for 
these changes in the porosity parameter. 
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RMS Error in Matching Cp on Walls 
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Figure 9. Velocity perturbation profiles on the top and bottom walls predicted by 
Navier-Stokes solver using a slip wall boundary condition. 



Figure 10. Root-mean-square errors as a function of the number of data points used to 
develop the method of singularities solution. Errors are normalized by peak pressure 
measured. 
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Figure 11a. Method of Singularities matches to pressure profiles on the wind tunnel 
walls for airfoil tests at a = 5 degrees. Uco = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 



Figure lib. Comparison of experimental and CFD pressure profiles on the airfoil for 
a = 5 degrees. U« = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 
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Figure 12a. Method of Singularities matches to pressure profiles on the wind tunnel 
walls for airfoil tests at a = 6 degrees. Uoo = 24. m/s. Re = 3.2 x 10 5 . h/c = 2.25. 
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Figure 12b. Comparison of experimental and CFD pressure profiles on the airfoil for 
a - 6 degrees. Uoo = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 
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Figure 13a. Method of Singularities matches to pressure profiles on the wind tunnel 
walls for airfoil tests at a = 7 degrees. Uo» = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 
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Figure 13b. Comparison of experimental and CFD pressure profiles on the airfoil for 
a = 7 degrees. Uoo = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 


51 






♦ Solid Wal Wind Tunnel Test Data 

CFD Solution using SokJ Wall Pressure Boundary 
Condition from the Method o< Singularities 
O Porous Wall Wind Tunnel Data 

CFD Solution using Porous Wall Pressure Boundary 
Condition trom the Method of Singularities 
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Figure 14b. Comparison of experimental and CFD pressure profiles on the airfoil for 
a = 8 degrees. Uoo = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 
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Figure 15a. Method of Singularities matches to pressure profiles on the wind tunnel 
walls for airfoil tests at a = 9 degrees. Uoo = 24. m/s, Re = 3.2 x 10 5 , h/c = 2.25. 



Figure 15b. Comparison of experimental and CFD pressure profiles on the airfoil for 
a = 9 degrees. Uoo = 24. m/s, Re = 3.2 x 10^, h/c = 2.25. 
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Figure 16a. Method of Singularities matches to pressure profiles on the wind tunnel 
walls for airfoil tests at a = 10 degrees. Uoo = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 



Figure 16b. Comparison of experimental and CFD pressure profiles on the airfoil for 
a = 10 degrees. Uoo = 24. m/s. Re = 3.2 x 10 5 , h/c = 2.25. 
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Figure 17b. Comparison of experimental and CFD pressure profiles on the airfoil for 
a = 1 1 degrees. Uoo = 24. m/s, Re = 3.2 x 10 5 , h/c = 2.25. 
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RMS Error in Matching Cp on Walls 
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Figure 18. Root-mean-square errors as a function of porosity parameter specified in 
method of singularities solution. Experimental data from a = 8 degrees, porous wall 
wind tunnel test 
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Figure 19. Method of singularities fits to experimental data using a porosity parameter 
of P = 0. and P = 1. Experimental data from the porous wall wind tunnel test at a = 8 
degrees. 
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Figure 20. Effect of using boundary conditions based on porosity parameters of P = 0. 
and P = 1. on computed airfoil pressure profiles for a = 8 degrees. Boundary conditions 
used are shown in Figure 19. 



Figure 21. Effect on singularity strengths of changes in porosity parameter specified in 
method of singularities. Method of singularities solutions were based on a = 8 degrees, 
porous wall wind tunnel test. 
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RMS Error In Matching Cp on Walts 
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Figure 22. Effect of changes in the position of singularities on the root-mean-square 
errors in matching pressure profiles for the a = 8 degrees, solid wall wind tunnel test 




Figure 23. Grid used in the a = 8 degrees simulation using viscous wall boundary 
conditions. 
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Figure 25. Pressure profiles on the airfoil from CFD calculations using specified 
pressure, no slip wall, viscous wall and free-air boundary conditions. 



Figure 26. Refined grid for a = 8 degrees test case. Grid dimensions are 350 x 105. 
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Figure 27. Comparison of results from 250 x 70 grid and 350 x 105 grid. Results are 
for solid wall test at a = 8 degrees. 
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Figure 28a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for solid wall test at a = 5 degrees. 



Figure 28b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 28a. 
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Figure 29a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for solid wall test at a = 6 degrees. 
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Figure 29b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 29a. 
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Figure 30a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for solid wall test at a = 7 degrees. 



Figure 30b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 30a. 


64 





Figure 31a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for solid wall test at a = 8 degrees. 



Figure 31b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 31a. 
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Figure 32a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for solid wall test at a = 9 degrees. 



Figure 32b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 32a. 
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Figure 33b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 33a. 
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Figure 34a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for solid wall test at a = 1 1 degrees. 



Figure 34b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 34a. 
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Figure 35a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for porous wall test at a = 5 degrees. 



Figure 35b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 35a. 
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Figure 36a. Comparison of method of singularities Fit to complete set of data and sparse 
data on the walls for porous wall test at a = 6 degrees. 
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Figure 36b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 36a. 


70 






Figure 37 a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for porous wall test at a = 7 degrees. 



Figure 37b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 37a. 
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Figure 38a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for porous wall test at a = 8 degrees. 



Figure 38b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 38a. 
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Fl gure 39b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 39a. J 
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Figure 40b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 40a. 
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Figure 41a. Comparison of method of singularities fit to complete set of data and sparse 
data on the walls for porous wall test at a= 11 degrees. 



Figure 41b. Comparison of CFD airfoil pressure profiles obtained using the boundary 
conditions developed from sparse and fine data, shown in Figure 41a. 


75 




Bottom Wall, 5-Hole Probe Data 
Top Wall, 5-Hole Probe Data 
Bottom Wall, CFD Result 
Top Wall, CFD Result 


Figure 42. Normal velocities near the walls, as measured by five-hole probe and 
computed by CFD code. Porous wall tests, a = 5 degrees. 
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Figure 43. Normal velocities near the walls, as measured by five-hole probe and 
computed by CFD code. Porous wall tests, a = 8 degrees. 








Figure 44a. Upper wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 

singularities solution. U„ = 24. m/s, Re = 2.4 x 1()5, h/c = 3.0, a = 20 degrees Solid 
wall boundary condition. ucgrees. solid 
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Figure 44b. Upper wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 

singularities solution. U« = 24. m/s. Re = 2.4 x 10 5 , h/c = 3.0, a = 20 degrees Solid 
wall boundary condition. 
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Figure 45a. Lower wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees Solid 
wall boundary condition. 
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Figure 45 b . Lower wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 
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Figure 46. Side wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. U©o = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 
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Figure 47a. Upper wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. Uo© = 24. m/s, Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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Figure 47b. Upper wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10 5 , h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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Figure 48a. Lower wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. U«> = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 








Figure 48b. Lower wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. Uoo = 24. m/s, Re = 2.4 x 10 5 , h/c = 3.0, a = 20 degrees. Porous 
wall boundaiy condition. 
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Figure 49, Side wall experimental data and method of singularities fit to the 
experimental data, using 16 singularities and 342 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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x/c 

Figure 50a-c, Locations of the complete set of 342 measurements on the top, bottom 
and side walls. 
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Figure 51a-c. Locations of the reduced set of 223 measurements on the top, bottom and 
side walls. Measurements are distributed along the centerlines of all three walls and 3 
cross-planes. 
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Figure 52a-c. Locations of the reduced set of 136 measurements on the top, bottom and 
side walls. Measurements are distributed along the centerlines of all three walls and 1 
cross-plane. 
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Figure 53a-c. Locations of the reduced set of 38 measurements on the top, bottom and 
side walls. Measurements are distributed along the centerlines of all three walls and 1 
cross-plane. 
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Number of Points Used to Determine Singularity Strengths 

Figure 54. Effect on root-mean-square errors of using various numbers of singularities 
and various numbers of data points to develop method of singularities solutions. Results 
are shown for solid wall test at a = 20 degrees. 



Number of Points Used to Determine Singularity Strengths 

Figure 55. Effect on root-mean-square errors of using various numbers of singularities 
and various numbers of data points to develop method of singularities solutions. Results 
are shown for porous wall test at a = 20 degrees. 
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Figure 56a. Upper wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. U«> = 24. m/s. Re = 2.4 x 10 5 , h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 
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Figure 56b. Upper wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. U©o = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 
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Figure 57a. Lower wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. U<» = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 
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Figure 57b. Lower wall expenmental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 












Figure 58. Side wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Solid 
wall boundary condition. 
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Figure 59a. Upper wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10 5 , h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 


97 




y/W = .9167 


y/W = .8056 


• 8 


• i • « 



• • 


• • 



y/W = .6944 


y/W = .5833 





• • 8 * 


y/W = .5278 


• Experimental Data 
O Method of Singularities Fit to Data 


Figure 59b. Upper wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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Figure 60a. Lower wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uoo = 24. m/s, Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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Figure 60b. Lower wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uco = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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Figure 61. Side wall experimental data and method of singularities fit to the 
experimental data, using 7 singularities and 38 measurements to determine method of 
singularities solution. Uoo = 24. m/s. Re = 2.4 x 10^, h/c = 3.0, a = 20 degrees. Porous 
wall boundary condition. 
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RMS Error in Matching Cp on Walls 
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Figure 62 Root-mean-square errors as a function of the porosity specified in the method 
of singularities. Data shown is for a = 20 degrees, porous wall test. 
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Chapter 6 
Conclusion 


A method for describing porous wall boundary conditions based on sparse, 
nonintrusive measurements of flow quantities at the wall boundary has been developed. 
This method utilized a potential flow solution based on least squares matching of 
singularity strengths to measured experimental data. The consistency of imposing a 
pressure boundary condition based on this method with other means of describing a wall 
boundary condition has been demonstrated in computational studies. This method has 
been shown to provide a good description of the entire wall boundary condition even 
when as few as 8 measurements were used to develop the method of singularities 
solution. The ability of this method to predict the effects of changing porosity in a wind 
tunnel test has been demonstrated by simulation of experiments performed in a variable 
porosity test section. 

This method has been extended to three-dimensional porous wall testing. In 
three dimensions, the method of singularities allowed the porous wall boundary 
condition to be modelled without the need for normal velocity perturbation 
measurements and without the need for extensive calibrations of the wall. An 
experimental study has shown that the method was capable of capturing the trends in 
pressure profiles existing on the walls in three-dimensional porous and solid wall tests. 

In low speed tests, the method of singularities has been found to be rather 
insensitive to the value of the porosity parameter specified. This allowed for the 
porosity parameter to be found by means of least squares matching. While changes in 
the porosity parameter did have a strong effect on the singularity strengths, the effect on 
the overall match to the pressure profile was minimal. Since the porosity parameter was 
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not specified as a boundary condition in the CFD code, and no corrections were made to 
the wind tunnel data, this method eliminated many of the concerns associated with the 
classical linear porous wall boundary condition. This method also allowed for, a 
simulation of the entire flow field and direct comparison of the flow field to wind tunnel 
data without the need for corrections to the experimental data. 
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